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Abstract 

We perform Euclidean strong coupling expansions for Yang Mills theory on 
the lattice at finite temperature. After setting up the formalism for general 
SU(N), we compute the first few terms of the series for the free energy density 
and the lowest screening mass in the case of SU(2). To next-to-leading order the 
free energy series agrees with that of an ideal gas of glueballs. This demonstrates 
that in the confined phase the quasi-particles indeed correspond to the T = 
hadron excitations, as commonly assumed in hadron resonance gas models. Our 
result also fixes the lower integration constant for Monte Carlo calculations of the 
thermodynamic pressure via the integral method. In accord with Monte Carlo 
results, we find screening masses to be nearly temperature independent in the 
confined phase. This and the exponential smallness of the pressure can be under- 
stood as genuine strong coupling effects. Finally, we analyse Pade approximants 
to estimate the critical couplings of the phase transition, which for our short 
series are only ~ 25% accurate. However, up to these couplings the equation of 
state agrees quantitatively with numerical results on A^ = 1 — 4 lattices. 



1 Introduction 



The study of QCD at finite temperatures and densities is of growing phenomenological 
interest for current and future heavy ion colhsion experiments, as well as for many 
astrophysical problems. Due to the interaction strength, the usual perturbative treat- 
ment by a series expansion in a small coupling constant fails for QCD. Perturbative 
failure persists in the quark gluon plasma phase, where the weak coupling series is 
only defined up to a certain maximal order (depending on the observable) before the 
well-known Linde infrared problem sets in [T]. The best one can do is to calculate an 
effective theory for the infrared modes with perturbatively calculable coefficients, such 
as dimensional reduction [2] or hard thermal loops [3J, and solve the effective theory 
for the soft modes on the lattice. The fully non-perturbative alternative are of course 
lattice Monte Carlo simulations. For a recent review of lattice results, see [1] . 

On the other hand, for temperatures below the plasma transition, no analytic ap- 
proaches starting from the QCD Lagrangian are available at all. Even lattice sim- 
ulations of the equation of state are difficult in this regime, due to the exponential 
suppression of the pressure with hadronic masses at low temperatures. Moreover, nu- 
merical determinations of the pressure via the integral method [5] require to supply 
a lower integration constant, corresponding to the pressure at some low temperature. 
Being unknown from first principles, this constant is usually set to zero by hand, based 
on its assumed exponential smallness. A successful description of lattice data below Tc 
is given by the hadron resonance gas model [6j, which requires expansive experimental 
knowledge of the hadron spectrum as well as some modelling. 

In the present work we fill these gaps with an analytic treatment of the low temper- 
ature phase by means of a strong coupling expansion of the lattice pure gauge theory. 
Contrary to weak coupling expansions, strong coupling expansions are known to be 
convergent series with a well-defined radius of convergence. In the early days of lattice 
gauge theory they were used to get analytical results for some physical quantities of 
interest, such as glueball masses [7] or the energy density of lattice Yang-Mills theories 
[8l[9]. These calculations were done at zero temperature, i.e. for lattices with infinite 
spatial volume and temporal extent Nf. To our knowledge, strong coupling expan- 
sions of the thermal partition function have not been considered beyond the infinite 
coupling limit, thus neglecting gauge fiuctuations. In the Euclidean framework, an 
effective theory for Wilson lines can be constructed in the strong coupling limit [TU] . 
to be analysed by mean field methods. The only series we are aware of is for the 
temperature-dependent string tension [11]. In Hamiltonian approaches, an effective 
Hamiltonian for the strong coupling limit is constructed, and the partition function 



1 



still has to be solved for by some other means. For the pure gauge theory, a deconfine- 
ment transition was predicted in this way [12]. More recent applications are to systems 
with fermions at finite density, e.g. [13]. For a review of early work and references, see 

In this work, we calculate Euclidean strong coupling series for the free energy density 
and screening masses in SU(2) pure gauge theory with an infinite spatial volume and 
a compactified temporal lattice extent Nt. In this way we can analytically study finite 
temperature effects in the confined phase. The physical deconfinement phase transition 
then corresponds to a finite convergence radius of the series, which we try to estimate 
from the behaviour of the coefficients. In Sec. [2] we set up the formalism of computing 
strong coupling series for the free energy at finite temperature. Sec. E] gives the explicit 
series for the SU(2) gauge theory and discusses how to leading order it coincides with 
that of a non-interacting glueball gas, as well as estimates for the radius of convergence 
and the phase transition. We also compare our results to Monte Carlo simulations. 
Sec. m discusses strong coupling series of screening masses before we conclude in Sec. O 



2 The free energy density as strong coupling series 



Consider SU(N) Yang-Mills theory on an f2 = A^^ x A^^ lattice with lattice spacing a 
and the Wilson action. Its partition function is given by 




Here, is the coupling constant of the corresponding continuum field theory, and the 
elementary plaquettes are given in terms of link variables as Up{x) = U^{x)Uy{x + 
a{l)Uj^{x + aP)Ul{x). The finite space-time box corresponds to a physical volume 
V = {aNs)^ and a temperature T = l/{aNt). We are interested in the thermodynamic 
limit and always consider the situation A^^ ^ oo. 

To our knowledge, Euclidean strong coupling expansions have so far only been ap- 
plied to the T = or Nt ^ oo situation. Formally, the partition function Eq. ([T]) 
still resembles a thermal system in this limit, with the lattice gauge coupling /3 playing 
the role of inverse temperature. Correspondingly, strong coupling expansions in small 
/3 are often termed 'high temperature expansions' in the literature [15], even though 
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the physical temperature of the system is zero. We shall not repeat the derivation 
of the strong coupling expansion here, but only give some central formulae required 
to fix the notation. For more details, see [15] and references therein. To start with, 
the exponential of the action is expanded in group characters Xr{U) = Tr Dr{U), with 
Dr{U) a specific irreducible representation matrix of U with dimension dr, 



The coefficients of the character expansion Cr{l3), and hence the effective expansion 
parameters ar(/5), can be expanded in powers of /5 to yield the desired series in g~'^. 
The series can then be reorganised as a sum of graphs G with contributions ^{G), 

Z = 4''J2^{G), $(G)= f DUl[d,^ar,XrAU) = 1[^{X.^^ (3) 

G peG i 

which factorise into disconnected components Xj, called polymers. Finally, using the 
formalism of moments and cumulants, one arrives at a cluster expansion for our quan- 
tity of interest, the free energy density, 

/^-ilnZ = -61nco(/3)-i J] a{G)l[<!>{X,r . (4) 

The sum is over all clusters C, which are defined as connected polymers Xi, and Ui 
denotes the multiplicity of a particular polymer in a cluster. The combinatorial factor 
a(C) is given as 

/r^\ [Xi, ■ ■ ■ , Xi, X2, . . . , X2, . ■ . , Xk] , . 

(^{Gj = j — j j (5j 

and equals 1 for clusters G which consist of only one polymer Xi. The so-called 
cumulant [ ] can be expressed in terms of moments < > 

[a, . . . , C] = - 1)! < a, . . . , /3 > ■ ■ ■ < 7, . . . , 5 > (6) 

p 

where n is the number of factors on the right hand side and the sum goes over all 
partitions P. The moments are defined in such a way, that 

{1, if every pair Xi^Xj is disconnected 
U, otherwise 

This implies the equivalence between non-zero cumulants and connectedness of graphs, 
[Xi, ...,Xn]^0 ^ Xi U ■ • • U X„ is connected. (8) 
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The contributing polymers Xj have to be objects with a closed surface since 



dUXriU) = 6r,0. (9) 

This means the group integration projects out the trivial representation at each link. 
Group integrals are calculated using integration formulae like 

J dUxr{UV)xs{WU-') = ^-^XriVW). (10) 

Note that, because of translation invariance and connectedness, the number of identical 
clusters at different positions is oc Q, so that the lattice volume drops out of Eq. (j4]) 
which thus has a finite thermodynamic limit. 

We now wish to apply this formalism to the case of non-zero physical temperature, 
which is realized by keeping A^^ = l/{aT) finite. The free energy density related to a 
physical temperature T is defined as 

- = --lnZ = ^\nZ = f(Nt), (11) 

i.e. it corresponds to the previously defined free energy density evaluated at finite 
Nt- Evidently, all formulae above remain unchanged in this case. The only effect of 
finite with periodic boundary conditions is to change the set of contributing graphs 
{Xi}. The physical free energy is then obtained by subtracting the formal free energy 
f{Nt = oo), which renormalises Eq. fill I) analogous to a subtraction of the divergent 
vacuum energy in the continuum, 

/(iV„/5) = /(iV„/5)-/(oo,/3). (12) 



2.1 Classification of graphs 

Because of the difference in Eq. f|T2l) . those graphs contributing in the same way to 
f{Nt) and f{oo) drop out of the physical free energy. This is true for all polymers 
with time extent less than Nf. The calculation thus reduces to graphs with a temporal 
size of Nt on the finite Nt lattice, and graphs spanning or extending A^^^ on the infinite 
lattice. Such graphs contribute either to f{Nt) or to f{oo) (and in some cases to both), 
and hence to the difference in Eq. [121 It is therefore clear from the outset that the 
strong coupling series for the physical free energy starts at a higher order than the 
formal zero temperature free energy. Moreover, the order of the leading contribution 
depends on N^. 
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Figure 1: Graphs that appear on the different lattices. Left: Leading order tube for 
Nt = 4. Middle: First graph vanishing in the case = 4. Right: Graph contributing 
on both, Nt = 4, CO, lattices. 



The lowest order graph existing due to the boundary condition on the finite Nt 
lattice, but not on the infinite lattice, is a tube of length A^^^ with a cross-section of one 
single plaquette, as shown in Fig. [1] (left). It forms a closed torus through the periodic 
boundary and thus gives a non-vanishing contribution, which is easily calculated to 
be $(Gi) = a^^*, where the subscript '/' indicates the fundamental representation. 
We need to sum up all such graphs on the lattice. There are three spatial directions 
for the cross section of the tube, giving a factor of 3. Translations in time take the 
graph into itself and do not give a new contribution, while we get V^{Gi) from all 
spatial translations. Together with the 1/Q in Eq. Hlthis gives a factor of 1/iVf. The 
contribution of all tubes with all plaquettes in the fundamental representation is thus 

HGi) = l-afN (13) 

which is - up to a sign - also the leading order result for the physical free energy. For 
SU(N) with > 3 we have an additional factor of 2 because there are also complex 
conjugate fundamental representations. 

On the other hand, the same tube with both ends closed off by additional plaquettes 
as in Fig. [1] (middle), contributes with $(G2) = '^/'^/^*^^ = oo, but not on a 

finite Nt lattice. This is because the boundary plaquettes get identified as one doubly 
occupied plaquette, which is not an allowed graph in the expansion. Therefore, $(6*2) 
counts with a negative sign relative to $(6*1) towards the physical free energy. Here, 
translations in time do produce a new graph, so the total contribution is 

$(G2) = -344^*+^ (14) 

Fig. [1] (right) shows a variation of these basic graphs, a cluster composed of two 
double cubes. This is an example of a graph spanning Nt which contributes to both 
the finite and infinite Nt lattices in the same way, thus cancelling out in the physical 
free energy. (Note that the corresponding tube obtained without the plaquettes at 
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Figure 2: Graphs contributing to the higher order terms of the series. 



Nt vanishes, it would correspond to a single polymer with doubly occupied plaquettes 
at the slit). However, for similar clusters composed of more than two polymers, this 
cancellation in general no longer holds because of different assignments of combinatoric 
factors a(C) in the two cases. 

2.2 Corrections to basic polymers 

For fixed Nt, starting from the basic leading order polymers discussed in the last section, 
one can now build up the corrections by adding decorations on each of them. These 
can be either geometric, by adding additional fundamental representation plaquettes 
as in Fig. [2] (left), by inserting plaquettes in higher representations as in Fig. [2] (middle) 
or by adding a whole new polymer as in Fig. [2] (right). Of course, these modifications 
can be combined. Adding plaquettes in higher representations is possible only if at 
each and every link the Clebsch-Gordan series of the representation matrices contains 
the trivial representation, due to Eq. Qj. 

We have already seen in Sec. 12. II that the order in a/, to which the graphs contribute, 
depends on A''^. Thus, the relative importance of different types of graphs changes with 
Nt. For our example Nt = 4 considered in Sec. 12.11 the leading correction to the basic 
polymer on the Nt lattice is the insertion of one higher representation plaquette inside 
the tube, contributing an additional factor ~ aj-, whereas the lowest order geometric 
decoration is a shifted plaquette adding a factor ~ a^. By contrast, on an Nt = I lattice 
contributions to the leading order correction ~ come from the basic polymer with a 
cross sectional perimeter of six links. In general, geometric decorations enter earlier the 
lower Nf is. Thus the summation of basic polymers and their decorations contains the 
complete result to some fixed order O only for lattices Nf > A*"^*^, with some N^ which 
is obviously growing with O. For lattices Nt < N^ there are additional geometric 
decorations contributing to 0{a^). In this work, we have calculated corrections to a^^* 
through O(a^) for which A^*^ = 5. 
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2.3 Summing basic polymers and their corrections for SU(2) 

The contribution of all graphs of length A^^ without geometric decorations can be 
summed up in closed form. To do this we note that an additional plaquette in a 
representation r as in Fig. [2] (middle) gives an additional factor of drttr. In the case 
of SU(2), for which we will describe the calculation, the only possibility is the j = 1 
representation, because of Eq. ([9]) and 

11 11 

-®-=0©l ^ -©-©1=0©1©1©2. (15) 

The expansion parameters of the lowest representations are given by modified Bessel 
functions and can be expanded in powers of the lattice coupling, 

" - !|i 4'' -95''' + 1^6'''^- SSO'^' + 

For the following it is convenient to introduce the combination c = 1 + 3f — and 
to use M = a/ as the effective expansion parameter instead of j3. It is well known from 
expansions at zero temperature that apparent convergence is better for the series in u 
[3, [9], and we observe the same phenomenon here. 

On the A''^ lattice we can have < k < Nt additional plaquettes at A^^ places which 
can be distributed in (^') ways. Summing over all possible distributions gives 

$1 = <i>(Gi) J2{k) = ^^^'^ + • 

k=0 ^ ^ 

We can also add slits to get graphs as in Fig. [T] (right), consisting of more than one 
polymer. Each slit gives a factor d'ju'^ = Av?. The minimum number of slits, i, is 2 
and the combinatorial factor for such graphs is 

a{C) = {^ - l){-ir\ (18) 

Summing these graphs with possible j = 1 plaquettes at the remaining places we get 

= *(GoE(T*)('-i)(-i'"'M"Er70'^"'' 

- (1 + ?>vf' + WAr,c^*-ij . (19) 
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Of course, we can make the same insertions to the graph Fig. [T] (middle) on the 
infinite lattice with the difference that we have at most Nf — l places to add plaquettes. 
The combinatorial factor now reads 

a{C) = (20) 
where i is again the number of slits and, in this case, it is unrestricted, giving 

= $(6*2)0^'-' 

= ^Gi){-Au^Nt)c^'-\ (21) 
For the final result we have to add the different pieces and get 

$ = $l + $2 + $3 = A^4Ar,^JV,_ (22) 

In higher gauge groups, the summation proceeds in a similar fashion with some slight 
modification due to the fact that there are also complex conjugate representations. 



3 The free energy density for SU(2) 

Let us now give our central results for the gauge group SU(2). For Nt > 5, the basic 
polymers and their decorations can be summed up to give the general result 



(23) 

For Nt = 1 — 4 there are additional geometric decorations, while some graphs contained 
in the previous result do not contribute on those short lattices. We find 

./X 4 fi 10913 s 968642 ^ ^ 12^ / 

f{l,u) = -3u'-16u'-^^u'-^^u'' + 0{u'') (24) 

fi2,u) = -^«« + 6«^°-55«- + ^«--^^^^«^^ + OK«) (25) 

f{3,u) = -u'^ + 6u''-50u''+^-^u''~^-^u^' + O{u^^) (26) 

/(4,m) = -^^16 + 6^18 _ 5g^2o ^ 51376 22 _ 2402453 ^4 q / 26^ ^^7) 
■'^ ' ' 4 135 810 \ I ^ ' 

The fact that only even powers of u appear is due to the reality of the SU(2) repre- 
sentations. Note that for Nt = 1, our series is shorter since we did not compute the 
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numerous geometric decorations at the next order. From these expressions all other 
thermodynamic quantities of interest can be constructed, in particular the pressure 
and energy density, respectively. 

Since the partition function is not directly measurable in Monte-Carlo simulations, the 
pressure is usually obtained by the integral method [5], where the expectation values 
of derivatives are computed and then integrated numerically, 

= Nt d/3'[6(Tr?7pO)-3(Trf/; + Tr?7;)], 

/3o J 1^0 

where (Trt/^ ) denotes the plaquette expectation value on symmetric (T = 0)lattices, 
Nt = Ns oo, and (Trf/^''') are those of space-time and space-space plaquettes for 
Nt < Ng. The lower integration limit is usually set to zero by hand, arguing with an 
exponentially small pressure in the low temperature regime. Our results justify this 
assumption from first principles and allow to fix that value if desired. 




3.1 The free energy density from an ideal glueball gas 

In weak coupling expansions of the pressure, the leading term is the well-known Stefan- 
Boltzmann limit, describing a non- interacting gas of the constituent particles. It is 
now interesting to ask how the QCD pressure can be interpreted in the strong coupling 
regime. From the Wilson action it is clear that the strong coupling limit is also non- 
interacting. However, as we have noted already, in this limit the pressure is zero. 
Considering strong but finite couplings, and recalling the first orders of the T = 
glueball mass calculations for SU(2) [7], 

o 2 98 4 20984 . 151496 . 

m(A++) = -4 In M + 2u^ u^ 

^ ^ ' 3 405 243 

.1 ^2 26 4 13036 28052 c 
m(E++) = -A\nu + 2u^ - —u' + ^^u' - ^^u\ (29) 

we observe that through second order in the expansion the free energy can be written 
as 

/(AT,, u) = --^ re-™(^^^)^* + 2e-"^(^^^)^* + 0{u^)] . (30) 

Here, the prefactors 1 and 2 before the exponentials correspond to the number of 
polarisations of the respective glueball states. Note that higher spin states start with 
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~ 6 Inw [TB], thus contributing to the order ~ u^^* or higher in the free energy. Hence, 
through two non-trivial orders our result is that of a free glueball gas, modified by 
higher order corrections. This is a rather remarkable result. It allows to see from a first 
principle calculation that the pressure is exponentially small in the confined phase, and 
that it is well approximated by an ideal gas of quasi-particles which correspond to the 
T = hadron excitations. While this result might be expected on phenomenological 
grounds, it is nice to see it demonstrated by an explicit calculation. This gives a 
quantum field theoretical explanation for the otherwise plausible success of the hadron- 
resonance-gas model in reproducing the confined phase equation of state [6]. 



3.2 Series analysis and phase transition 

Strong coupling or high temperature expansions have been worked out to high orders 
in many spin models, where various tools of series analysis can be applied to improve 
convergence or extract additional information from the behaviour of the coefficients P, 
[T7] . In this section we explore some of these possibilities with our series. In particular, 
we are interested in the radius of convergence of the strong coupling series. SU(N) 
pure gauge theories at finite temperature have a true order parameter for confinement, 
the Polyakov loop, and therefore there is a non-analytic phase transition separating 
the confined from the deconfined phase. The associated critical coupling jSc limits 
the radius of convergence of the strong coupling series, provided there are no other 
singularities Ps in the complex /3-plane with \j3s\ < Pc- 

For our later comparison with Monte Carlo results, it is particularly convenient to 
consider the energy density. The curves for consecutive orders in the strong coupling 
expansion for Nt = 2 are plotted in Fig. [3] (left). For (3>1 convergence rapidly becomes 
poor, announcing the proximity of the convergence radius. The most straightforward 
way to estimate the radius of convergence r„ from the nth order series would be via 
the ratio test of the series coefficients f2n, 



f: 



2n 



f: 



2n+2 



1/2 

(31) 



However, our series Eqs. fl23ti27p are still rather short and no convergence in r„ is 
visible. Moreover, a singularity on the real axis requires the coefficients asymptotically 
to come with equal signs. This is clearly not the case in our expressions for iVt > 1, 
which suggest a nearby imaginary singularity. 

A much better tool for our purposes is the analysis of Fade approximants to a 
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Figure 3: The energy density in the confined phase on Nf = 2. Left: subsequent orders 
of the strong couphng series. Right: Reconstructed from the L + M = 4, 6, 8 Fade 
approximants to the strong couphng series of Dc{u), Eqs. (l35l) . ( l38l) . 



function constructed from its series expansion. These are the rational functions 



ag + aiu + ■ ■ ■ + aiu^ 
l + biU-\ h b^ju^ 



\L,M](U) = ^ ^ ry, 



with coefficients Oj, 6, chosen such that they reproduce the power series of the function 
of interest to the degree L + M. As rational functions, Fade approximants are known 
to give good estimates of isolated pole singularities, whereas branch cuts or algebraic 
singularities are less well reproduced (for a detailed discussion, see [E]). Furthermore, 
they give access to several singularities, rather than just the nearest one. At the finite 
temperature phase transition the free energy f{Nt,Uc) with Uc = u{j3c) is continuous, 
with a discontinuous first or second derivative, depending on the order of the transition. 
This type of singularity is difficult to model for Fade approximants. Instead, the 'heat 
capacity' 
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C{N„u)=u'^J{N„u) (32) 

diverges at the phase transition as C{u) ~ {uc — m)" with a critical exponent charac- 
teristic of the transition. Its logarithmic derivative 

Dc{Nt, u) = ^\n C{Nt, u) ^ (33) 

au Uc — u 

has a simple pole with residue a and is therefore best suited for an analysis by Fade 
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approximants. We thus consider the series 

1879249 





- -( 


1 + 


20 2 54791 4 

— u H u 

3 243 


Dci2,u) 


- -( 




45 2 6445 4 
28 196 


Dc{3,u) 


12 

u 




91 2 4573 4 
" 66^ + 242 " 




16 

u 




51 2 29791 

M H u 

40 1800 



486 



(34) 



150331 31831541863 o\ , , 

+ (35) 

92610 21003948 ) ^ ' 

2653298 « 114561591157 o\ 

H M*^ (36) 

59895 106732890 ) ^ ' 

1262057 1055884297 o\ 

+ (37) 

27000 1215000 J ^ ' 

and model the full functions by Pade approximants. Using the formulae 

C{u) = exp j duDciu), ^f{u) = j du^, (38) 

the energy density can be reconstructed from the Pade approximants to Dc{u). As an 
example, we show the L + M = 4, 6, 8 approximants for the Nt = 2 lattice in Fig. [3] 
(right). Clearly, this sequence of approximants shows improved convergence compared 
to that of the bare series. 

The singularities in Dc{u), indicated by zeroes of the denominator, together with the 
zeroes of the resulting approximants are shown in Table [T] for Nt = 2, 4, respectively. 
Singularities in the immediate neighbourhood of a zero of the same approximant are 
typically artefacts and unstable under variation of the approximant. However, for Nt = 
2 several approximants show poles around (3 = 1.5 without zeroes in the immediate 
vicinity, indicating that the full function indeed has a singluarity on the real axis in 
this region. For A^^^ = 4, on the other hand, the pole near /? = 1.5 is accompanied by 
a zero and not to be taken seriously. The next nearest pole on the real axis is instead 
around (3 = 2. A priori it is not possible to judge which approximants are better than 
others. The scatter in the results is thus a measure for the systematic error associated 
with the Pade procedure. Moreover, there is also a scatter between approximants based 
on different orders of the underlying series. With increasing order, the approximants 
display more and more singularities which should eventually accumulate near the true 
singularity structure. 

To take these systematic effects into account, we estimate (3c by averaging over 
the lowest lying real singularities obtained from the two highest order approximants, 
i.e. L + M = 4, 6 for A'f = 1 and L + M = 6, 8 for A"* = 2 - 4. To quantify the scatter 
due to the systematic uncertainties we quote (/5™"^' — /?™")/2 as an error estimate. 
The same procedure is followed for the residues, and the results are collected in Table 
[21 Since the series are still short, the predictions for the critical coupling are not yet 
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Pade 


Singularities 


Zeroes 


Residues 


[0,6] 


±1.8356 




0.07486 




±(0.2112 ± 1.7192z) 






[2,4] 


±1.6995 


±2.8474 


0.05656 




± 1.5674 z 






[4,2] 


± 9.9475 i 


±(1.1253 ± 1.2410i) 




[0,8] 


±1.6930 




0.06080 




± 1.4349 z 








±(0.9152 ± 1.6609i) 






[2,6] 


± 1.4893 


±1.7142 


0.02888 




± 1.5209 z 








±4.2104 






[4,4] 


±1.5502 


±1.9430 


0.03792 




±1.4997z 


± 2.2802 i 




[6,2] 


± 0.1308 i 


± 0.1308 2 








± (1.1278 ± 1.23842) 




Pade 


Singularities 


Zeroes 


Residues 


[0,6] 


±3.1022 




0.13287 




± 1.6250 i 






[2,4] 


±3.0636 


±3.0406 






± 1.6228 z 






[4,2] 


±1.9628 2 


±(0.9730 ± 1.7187i) 




[0,8] 


±2.0270 




0.05936 




± 1.4583 i 








±(1.0967 ± 1.5608i) 






[2,6] 


±2.8687 


± 0.4335 i 


0.11573 




±0.4335ii 








± 1.6505i 






[4,4] 


±1.5084 


±1.5762 






± 1.4023 i 


± 1.5976 i 




[6,2] 


± 0.8961 i 


± 0.8989 i 








±(1.2985 ± I.59I62) 





Table 1: Singularities, zeroes and residues of L ± M = 6, 8 Pade approximants for 
Nt = 2 (top) and Nt = 4 (bottom). Residues are only given for those singularities that 
enter the estimates for Pc, Table [2l 
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Nt 


a 




Pc (Monte Carlo) 


1 


0.061(38) 


0.92(15) 


0.85997(10) [18] 


2 


0.052(19) 


1.65(35) 


1.880(3) [19] 


3 


0.078(50) 


2.26(63) 


2.177(3) [19J 


4 


0.102(37) 


2.66(54) 


2.299(6) [T9| 



Table 2: Estimates for the critical coupling /?c and the critical exponent of the decon- 
finement phase transition. The exponent for 3d Ising universality is a = 0.12. 

very accurate, and those for the critical exponent even less so. Note, however, that a 
first order phase transition has a = 0, whereas a second order transition in the 3d Ising 
universality class has a = 0.12. Our results clearly favour the latter, especially as the 
lattice becomes finer. 

3.3 Comparison with Monte Carlo data 

It is now interesting to compare the results from the strong coupling series with Monte 
Carlo simulations. The thermodynamic quantity most easily accessible by Monte Carlo 
is the energy density, which is simply the expectation value of the plaquette, 

ei/3) = l^fiP) = (Trf/,)^, - {TtU,)^,=^, (39) 

where again the zero temperature (infinite Nt) piece is subtracted for renormalisation. 
As we have seen, in the low beta region of the deconfined phase, the corresponding 
values are exponentially small, and very high statistics runs are necessary in order to 
get significant results for a quantitative comparison. For the infinite volume vacuum 
lattice we have taken 12^^, and Ns = 12, Nt = 1, 2, 3, 4 for the finite T lattices. On the 
Nt = 2 lattice up to 1.5 x 10^ field configurations were generated to achieve sufficient 
accuracy, for the larger iV^'s this gets scaled down accordingly. 

We compare these data with the best estimate based on the strong coupling series, 
i.e. the Fade approximants to the highest available order in the logarithmic derivative 
of the heat capacity. Detailed results for Nt = 2,3 are shown in Fig. |H The different 
curves correspond to different approximants to the same order, and thus serve as an 
error band quantifying the uncertainties associated with the Fade procedure, thus 
giving a valuable error estimate. We observe quantitative agreement with the lattice 
data all the way up to the lowest estimates of Pc- For Nt = 1,4: we have checked at 
a few points that a similar picture obtains. Thus, the error estimate based on our 
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Figure 4: Comparison of Monte Carlo data for Nt = 2 (left) and Nt = 3 (right) with 
Fade approximants from the strong coupling series. 

Fade analysis appears to be reliable and announces the breakdown of the validity of 
the series. 



4 Screening masses 

Screening masses are defined by the exponential decay of the spatial correlation func- 
tion of suitable operators. An overview regarding definition, quantum numbers and 
numerical results can be found in [20]. Here we consider the colour-electric field cor- 
relator ( Tr FQj(x) Tr FQj(y) ), which is in the J^S"*" = 0^"*" channel (T denotes reflection 
in Euclidean time) containing the ground state and the mass gap. On the lattice, this 
corresponds to a correlation of temporal plaquettes, and the quantum numbers under 
the point group are A'^^. 

Temporarily assigning separate gauge couplings to all plaquettes, the correlator can 
be defined as jTj 



■In Z{P,P,P2: 



(40) 



/3i,2=/3 



At zero temperature the exponential decay is the same as for correlations in the time 
direction, and thus determined by the glueball masses, the lowest of which may be 
extracted as ^ 

m = - lim - In C{z). (41) 



The leading order graphs for the strong coupling series at zero temperature are shown 
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Figure 5: Graphs contributing to the lowest order of the expansion of the screening 
mass at vanishing and finite temperature. The correlated plaquettes are black. 

in Fig. [5] (left). This leads to the lowest order contribution: 

C{z) = Au^' = Ae-""''. (42) 

Thus, to leading order for the glueball mass is = —4 In u{l3). 

Now we switch on a physical temperature, i.e. keep the lattice volume compact in 
the time direction. As in the case of the free energy, we are here only interested in the 
temperature effects, i.e. in the mass difference 

Am(T) = m(T) - m(0) = - lim - [In C(T; z) - In C(0; z)] (43) 

with AC{T;z) = C{T;z) — C{0;z). A typical graph contributing in lowest order to 
this difference is shown in Fig. O (right). Summing up all leading and next-to- leading 
order graphs gives 

Am(r) = --Nt M^^*-^ c^* (1 + Au^), (45) 
3 

i.e. the screening masses decrease compared to their T = values. As in the case of 
the free energy, due to the difference only A/'rdependent higher orders contribute to the 
temperature dependence of screening masses. Again, the leading order result is generic 
for all SU(N) and quantum number channels. We conclude that in the confinement 
phase the lowest screening masses in each quantum number channel should be close 
to the corresponding zero temperature particle masses, with a significant temperature 
dependence showing up only near T^. This explains the findings of numerical investi- 
gations of the lowest screening mass in SU(3) gauge theory, which for temperatures as 
high as T = 0.97Tc see very little temperature dependence, Am(T)/m(0) > 0.83 [2T] . 
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5 Conclusions 



We have explored the possibihties of Euchdean strong couphng expansions for finite 
temperature lattice Yang Mills theories in the confined phase. The general formalism 
applying to zero temperature calculations can be taken over to this case, with the 
compact time dimension affecting the type and number of graphs that contribute to a 
certain quantity. As a consequence, temperature effects on the free energy density and 
screening masses appear only at an A^^-dependent higher order and vanish exponentially 
a,s f3 — 2N/g^ — > 0. This explains the numerically observed exponential smallness of 
the pressure and the near temperature independence of screening masses in the confined 
phase typical strong coupling phenomenon. 

We have explicitly calculated the first five terms of the series for the free energy 
density in the case of SU(2). To the leading two orders, the result agrees with that of 
an ideal glueball gas. This demonstrates that in the confined phase the quasi-particles 
indeed correspond to the T — hadron excitations, as is typically assumed in hadron 
resonance gas models. For lattices with A^ = 1 — 4, we have analysed the series with the 
help of Fade approximants and estimated the critical couplings Pc of the deconfinement 
phase transition. Since the series are still relatively short, those rcsTilts are not very 
accurate yet. However, within the estimated errors they are consistent with the critical 
couplings observed in Monte Carlo simulations. Moreover, up to the lowest estimated 
values of f3c the Fade approximants give a quantitative description of the Monte Carlo 
data for the equation of state. 

In conclusion, the strong couphng expansion offers valuable qualitative insight into 
temperature effects of the fully interacting, non-perturbative theory in the confinement 
phase. We are currently extending this work to the gauge group SU(3). 
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